f=[2;1];
A=[-4,-3;1,2];
b=[-6;3];
Aeq=[3,1];
beq=[3];
Ib=zeros(2,1);
[x,fval]=linprog(f,A,b,Aeq,beq,Ib)